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Abstract 

We investigate energetic stability and dissociation dynamics of water adsorption at the LaAlOa surface of 
the n-type LaAlOs/SrTiOa (LAO/STO) interface and its effect on electronic properties of the interface by 
carrying out first-principles electronic structure calculations. In an ambient atmosphere at room temperature 
the configuration of 1 monolayer (ML) of water molecules including 3/4 ML of dissociated water molecules 
adsorbed at the surface is found to be most stable, the configuration of 1 ML of dissociated water molecules 
is metastable. Water molecule dissociation induces a shift-up of the valence band maximum (VBM) of the 
LAO surface, reducing the gap between the VBM of the LAO surface and the conduction band minimum 
of the STO. For the LAO/STO interface with three LAO unit-cell layers, once the coverage of dissociated 
water molecules reaches 1/2 ML the gap is closed, the interface becomes metallic and the carrier density at 
the LAO/STO interface increases with increasing the coverage of dissociated water molecules. Our findings 
suggest two ways to control the conductivity at the LAO/STO interface: (I) insulator-metal transition by 
adsorbing a mount of water at the bare surface; (II) carrier density change by the transition between the 
most stable and the metastable adsorption configurations for 1 ML coverage in an ambient atmosphere at 
room temperature. 
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I. INTRODUCTION 



Two-dimensional electron gas (2DEG) at interfaces has been one of the most intriguing subjects 
of electronic material research. Because of its versatile and remarkable properties many attentions 
have been drawn from both fundamental researches and application studies. One of the key find- 
ings in this subject is that a 2DEG can be generated at the interface between two insulating metal 
oxides, LaAlOs/SrTiOa (LAO/STO)^ At the n-type LAO/STO interface between (LaO)+ and 
(Ti02)° layers, a high density (a few 10^^ cm"^) and high mobility 10^ cm^V~^s~^) electron 
gas has been realized at low temperature^. In addition, ferromagnetism^ and superconductivity^ 
were also reported at the same interface. As LAO and STO are two band insulators, finding the 
electron gas at their interface immediately triggered extensive researches to reveal the formation 
mechanism of the electron gas. Up to now, two mechanisms have been proposed and well ac- 
cepted. One is the polar catastrophe, in which the polar electric field in the LAO film shifts up 
the valence band maximum (VBM) of the LAO film layer by layer and finally leads to a charge 
transfer from the valence band of the LAO surface layer to the conduction band of the STO"^. It 
was found in Thiel et aUs experiment and confirmed theoretically that there is a critical thickness 
of the LAO film (4 unit cell) below which the interface is insulting and above which the interface 
becomes metallic^"-. The other mechanism is the role of oxygen vacancies, which were explic- 
itly demonstrated in experiments to exist at the STO side for the samples grown at low oxygen 
pressures^. Theoretical study showed that oxygen vacancies also readily form at the LAO sur- 
face and produce a carrier charge with a density of no more than half electron per u.c. area at the 
interface^. 

Recently, Cen et al reported a novel method to generate an electron gas at the LAO/STO 
interface^*i^. In their experiment, nanoscale control of so-called insulator-to-metal transition 
(IMT) at the LAO/STO interface was realized by an atomic-force-microscopy (AFM) tip scan- 
ning the LAO surface. A nanoscale conducting line can be formed at the interface by the AFM 
tip with a positive bias voltage writing on the LAO surface and the conducting line can be erased 
by the AFM tip with a negative bias voltage cutting the line. It is noticeable that this IMT only 
happens to the interface with three unit-cell (u.c.) layers of LAO. Subsequent experimental studies 
suggested that this IMT may be related to water adsorption at the LAO surface^»i^. To interpret this 
phenomenon several assumptions have been proposed. Bi et al assumed that the adsorbed water 
are dissociative and the AFM tip with a positive bias voltage removes the hydroxyls (OH~) from 
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the surface, as a result the left hydrogen atoms dope electron carriers at the interface However, 
first-principles electronic structure calculations have showed that regardless of the LAO thickness 
hydrogen adsorption on the LAO surface always induces a metallic interface^^, implying that the 
change at the LAO surface induced by the AFM tip scanning is not like Bi et aUs scenario. Xie et 
al speculated that water molecules are adsorbed at the LAO surface by a dipole-dipole interaction 
(the dipole of water molecule and dipole of the LAO film)^^, i.e., two hydrogen atoms of a water 
molecule approach to the LAO surface and the oxygen atom points to the outside. This will lead 
to an enhanced polarity of LAO film with adsorbed water compared to the bare LAO film, prob- 
ably resulting in a charge transfer from the surface to the interface. However, many studies have 
shown that water adsorption on oxide surfaces is via the formation of a chemical bond between 
the oxygen atom of the water molecule and the metal atom at the oxide surface^"— instead of 
a dipole-dipole interaction, basically ruling out Xie et aUs scenario. This is because that the en- 
ergy gain of the chemical adsorption on oxide surfaces, 1 ^3 eV per water molecule^»^»2^»21^ is far 
larger than that of the physical adsorption, about 0.2 eV per molecule^^. Up to now, although this 
nanoscale control of transport property at the LAO/STO interface has presented a very promising 
application on nanoscale memory device, called 'atomic force pencil and eraser—, its physical 
nature is still a puzzle. 



In this paper we present a first-principles study of electronic and structural characteristics of 
water molecules adsorbed at the LAO surface of the n-type LAO/STO interface, including sur- 
face atomic structures, thermodynamic stabilities, dissociation processes of water molecules, and 
electronic properties. We find that in an ambient atmosphere at room temperature the adsorp- 
tion configuration of one monolayer (ML) of water molecules including 3/4 ML of dissociated 
water molecules is most stable and the configuration of 1 ML of dissociated water molecules is 
metastable. The dissociation of water molecules at the LAO surface can reduce and even close 
the gap between the VBM of the LAO surface and the conduction band minimum (CBM) of the 
STO, leading to a metallic interface. Remarkably, for the ground-state and metastable adsorption 
configurations electron carrier densities at the interface show a large difference. Based on our 
findings we make an attempt to reveal the mechanism of the mysterious IMT found by Cen et al. . 
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II. COMPUTATIONAL DETAILS 



We carried out first-principles electronic structure calculations within the frame of density func- 
tional theory as implemented in the Vienna ab initio simulation package (VASP^. A generalized 
gradient approximation (GGA)^ for exchange-correlation functional and the projector augmented 
wave method22i22 with a kinetic cutoff energy of 400 eV for the plane wave basis set were used. The 
n-type LAO/STO interface was modelled into a (2x2) supercell consisting of four STO(OOl) lay- 
ers, three LAO layers, and a vacuum region of about 14 A. Various coverages of water molecules 
were located at the AIO2 -terminated surface. Dipole correction was employed to correct the errors 
of electrostatic potential, atomic force and total energy caused by periodic boundary condition^. 
The in-plane lattice constant of the slab was constrained at the calculated equilibrium value of 
STO bulk (a = 3.94 A). A (4x4x1) k-point grid was used for calculations of atomic structure 
optimization and a (6x6x1) k-point grid for electronic structure calculations. All atomic coor- 
dinates were fully relaxed until the atomic forces were less than 0.02 eV/A, except for the atoms 
in bottom two STO layers, which were fixed in their bulk positions. A first-principles molecular 
dynamics simulation at T=300 K was performed to search the atomic structure of 2 ML of water 
molecules adsorbed at the surface. The nudged elastic band method was used to calculate reaction 
paths of water molecule dissociation at the surface. 

III. THERMODYNAMIC STABILITY OF WATER ADSORPTION 

Figure [I] displays geometries of various adsorption configurations at different coverages, in 
which water molecules are either molecular or dissociated. For the molecular adsorption, for 
instance c(l/4, 0), the O atom of the water molecule sits on the Al atom forming an Al-O bond, and 
two H atoms approach to two closely neighbouring O atoms at the AIO2 layer. For the dissociative 
adsorption, for instance c( 1/4, 1/4), two distinct hydroxyl groups are produced, one of which 
originates from the water molecule and sits on the Al atom, the other (Ol-Hl) lies in the AIO2 
layer. For the configurations of no more than 1 ML coverage water molecules are bound to the 
surface by valence bonds or ionic bonds. Once the coverage exceeds 1 ML, as shown in c(2, 1/2), 
the extra molecules are loosely adsorbed at the surface by hydrogen bonds. 

Since the IMT observed in Cen et al. 's experiments happened in a humid ambient atmosphere at 
room temperature, to directly address the experimental results we first investigated thermodynamic 
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FIG. 1. (Color online) Top views of adsorption configurations of water molecules adsorbed at the AIO2- 
terminated surface of the n-type LAO/STO interface. The notation c(0, 6) denotes the adsorption configu- 
rations, in which 9 is the total coverage of water molecules and S is the coverage of dissociated molecules 
in unit of ML. Purple balls represent Al atoms, red O atoms, and grey H atoms. In c(2, 1/2) green dot lines 
represent hydrogen bonds. 
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FIG. 2. (Color online) Surface free energies of various adsorption configurations at the LAO surface as 
functions of the relative chemical potential of water vapor, Afia^Q. The top axis shows the relative humidity 
of the ambient atmosphere at 300 K corresponding to the relative chemical potential on the bottom axis. 
The free energy of the bare surface is taken to be zero. 
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stabilities of various adsorption configurations in an ambient atmosphere at room temperature by 
comparing their surface free energies. We calculated the surface free energies by employing the 
first-principles thermodynamics approacb^l, in which the surface is assumed to be in equilibrium 
with a humid ambient atmosphere characterized by the relative chemical potential of water vapor. 
The most stable configuration at a given relative chemical potential minimizes the surface free 
energy approximately defined as 

7 = ^4^^^^ ~ ^bare " ^HsO ' [^H20(gas) + A/XH2O PH20)] } • 

Here, E'ad is the total energy of the interface with ria^o adsorbed water molecules per surface area 
^bare IS thc total cucrgy of the bare surface, £'H20(gas) is the total energy of a gaseous water 
molecule, A/XH2o(r, Pr^o) the relative chemical potential of the water vapor at the temperature 
T and partial pressure pu^o. The vibrational free energy and configuration entropy contributions 
are neglected in the definition of surface free energy. E'ad, £'bare. and £'H2 0(gas) were obtained from 
our first-principles calculations. A/XH2o(r,l^H2o) follows from the ideal-gas relation 

A/XH,o(r,pH2o) = A/XH,o(r,p^,o) + kBTln{^) . 
At standard pressure p^^Q=l bar the relative chemical potential can be derived using 

where the enthalpy and entropy differences are obtained from thermochemical tables^^. The am- 
bient atmosphere at room temperature sets an upper limit for A/XH20 = — 0.57 eV per water 
molecule, which corresponds to the relative chemical potential of water molecules in a water- 
saturated atmosphere (100% relative humidity) at room temperature. The chemical potential at 
any other relative humidity can be obtained from the above ideal-gas relation. Figure |2] display s 
the surface free energies of various adsorption configurations varying with the relative chemical 
potential or corresponding relative humidity at room temperature. In the range of relative hu- 
midity varying from 0.01% to 100%, which nearly covers all possible relative humidities of the 
ambient atmosphere, the most stable configuration is c(l, 3/4). The configurations with a lower 
water coverage, such as c( 1/4, 1/4), can only exist in an extremely dry environment. Noticeably, 
the configuration c(2, 1/2) with 2 ML of water has a higher free energy than c(l, 3/4) and even 
than c(l, 1/2). This implies that once the coverage reaches 1 ML the surface will not adsorb 
water molecules from the atmosphere anymore, instead, water molecules readily remain in the 
atmosphere. 
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IV. ADSORPTION OF 1/4 ML OF WATER MOLECULES: DISSOCIATION PROCESS AND 
ELECTRONIC PROPERTIES 
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FIG. 3. (Color online) (a) Dissociation process of 1/4 ML of water molecules adsorbed at the AIO2- 
terminated surface of the LAO/STO interface. Reaction coordinate is the distance doi-Hi between the 
hydrogen atom HI and the oxygen atom OL Insets are perspective views for molecular, intermediate and 
dissociative configurations, (b), (c) Projected densities of states (PDOSs) for Ol atom in (b) the molecular 
(mol) and (c) dissociative (dis) adsorption configurations of 1/4 ML coverage, respectively. The vertical 
lines denote the Fermi level. The VBM of the STO bottom layer is taken to zero (see Fig. |4]). The x and 
y axes are approximately along the direction of Al-Os-Al chain and that of Al-Ol-Al chain shown in the 
insets in (a). 



For simplicity and clarity we first explored the dissociation process of 1/4 ML of water 
molecules adsorbed at the LAO surface and its effect on electronic properties of the interface 
system. Figure Oa) displays a dissociation process without energetic barrier for 1/4 ML of water 
molecules at the LAO surface, in which the binding energy is defined as £'b = £'ad — ^bare — 
^HsO • ^H20(gas)- The intermediate state involves a four-member ring Al-Ow-Hl-Ol. Generally, 
a water molecule dissociation at metal oxide surfaces is associated with a breaking of the metal- 
oxygen bond at the surface^»^»2^»21i22^ implying that the dissociation occurs more easily at the 
surface with loose metal-oxygen bonds. Indeed, the Al-O bond at the LAO surface is far weaker 
than that of a-Al203. This fact is reflected by the bond lengths in these two materials, 1.95 A 
(in-plane) for the former and 1.70 A for the latter^^. This implies that the energy cost to stretch the 
Al-Ol bond in the four-member ring is relatively low and can be compensated by the formation 
of hydroxyl Ol-Hl bond. To reveal the driving force of water molecule dissociation at the LAO 
surface we checked on the change of electronic states of Ol atom in the molecular and dissociative 

7 



adsorptions. In the molecular adsorption [Fig. Hb)] the orbital of Ol is a non-bonding state 
lying near the top of valence band. While in the dissociative adsorption [Fig. Oc)] the bonding 
state of the hydroxyl Ol-Hl mixing part states of and Py orbitals lies at the bottom of valence 
band. The lower surface energy of the dissociative adsorption just arises from the falling down of 
band energies of these O p states. Therefore, we see in Fig. |2lthat for the configuration of 1 ML 
coverage the surface free energy decreases with increasing the coverage of dissociated molecules 
up to 3/4 ML. The energy rising of 1 ML of dissociated water molecules is due to a large structural 
distortion at the surface. 
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FIG. 4. (Color online) (a) Layer-resolved projected densities of states (LPDOSs) for the configurations 
of 1/4 ML coverage. Grey-filled lines represent the LPDOSs for the molecular adsorption, red solid lines 
represent the LPDOSs for the dissociative adsorption. The PDOSs of water molecule are amplified four 
times for clear vision, (b), (c) Projected densities of states (PDOSs) for Os atoms in (b) the molecular (mol) 
and (c) dissociative (dis) adsorption configurations of 1/4 ML coverage, respectively. The vertical lines 
denote the Fermi level. The VBM of the STO bottom layer is taken to be zero, (d), (e) Potential energies 
averaged along the y axis over the slab of 2 A centering at Os atom in (d) the molecular and (e) dissociative 
adsorption configurations of 1/4 ML coverage. 
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Figure lit a) shows LPDOSs of molecular and dissociative adsorptions of 1/4 ML coverage. We 
note that water molecule dissociation gives rise to a remarkable change of electronic states of the 
surface layer, compared to the molecular adsorption the VBM of the LAO surface layer shift up 
about 0.2 eV in the dissociative adsorption. PDOSs in Figs. |4tb) and Uc) show that this shift- 
up is mainly due to the shift-up of the orbital of the O atoms like Os atom near the hydroxyl 
Ow-H. We ascribe this shift-up to an extra repulsive Coulomb potential generated by the negative- 
charged hydroxyl Ow-H. As shown in Fig. ffid) andllte), compared to the molecular adsorption, 
in the dissociative adsorption the potential energy over Os atom rises up significantly, implying 
that the orbital of Os atom lies at a higher position. Moreover, our calculations show that as 
more water molecules dissociate the VBM of the LAO surface layer shifts up more. A coverage 
of dissociated water molecules as high as no less than 1/2 ML can even shift up the VBM of the 
LAO surface layer by about 0.6 eV, consequently closing the gap between the VBM of the LAO 
surface layer and the STO CBM and thus leading to a metallic interface. 



V. ADSORPTION OF 1 ML OF WATER MOLECULES: DISSOCIATION PROCESS AND ELEC- 
TRONIC PROPERTIES 
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FIG. 5. (Color online) Dissociation process of 1 ML of water molecules at the LAO surface. Reaction 
coordinates are the distances dni-Oi between the hydrogen atom Hz and the oxygen atom Oi, z = 1, 2, 3, 4. 

Dissociation processes of 1/2, 3/4 and 1 ML of water molecules adsorbed on the LAO sur- 
face were also investigated. We found that for all of them there is an up-limit coverage of water 
molecules that dissociate barrierlessly, 1/2 ML. Here we are most interested in the case of 1 ML 
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coverage because c(l, 3/4), a configuration of 1 ML coverage, is the thermodynamic ground state 
in an ambient atmosphere at room temperature and as shown below the change of electronic prop- 
erties of the LAO/STO interface during the dissociation process of 1 ML of water molecules can 
be compared with the experimental results. Figure |5] shows the energy profile of a reaction path 
for the dissociation of 1 ML of water molecules, in which four water molecules dissociate one by 
one. We note that the barrierless dissociation of preceding two water molecules induces a large 
distortion at the surface. This distortion largely increases the distance between the hydrogen atom 
H3 of the third water molecule and the oxygen atom 03 at the LAO surface, finally giving rise 
to an energy barrier of 0.1 eV in the dissociation path of the third water molecule. Furthermore, 
the dissociation of the fourth water molecule undergoes an even higher barrier of 0.41 eV and the 
association of this dissociated molecule undergoes a barrier of 0.17 eV. 



To clearly see the effect of the dissociation of 1 ML of water molecules on the electronic prop- 
erties of the LAO/STO interface we plot LPDOSs in Fig. |6] (a) for 1 ML coverage. PDOSs of 
water molecules and the LAO surface layer present apparent differences in the molecular and dis- 
sociative adsorptions. In the fully molecular adsorption three molecular orbitals of water molecule 
corresponding to three peaks ranging from -7 to -2 eV are nearly well preserved. This molecular 
adsorption barely affects the electronic states of the LAO surface, the band gap between the VBM 
of the LAO surface and the CBM of the STO is 0.54 eV, approximately same as the value of the in- 
terface with a bare surface, 0.56 eV. In contrast, water molecule dissociation significantly modifies 
the electronic states of the LAO surface and the interface. Accompanying with the dissociation 
the VBM of the LAO surface layer shifts up considerably and a surface state originated from the 
hydroxyls on Al atoms appears at the VBM. Strikingly, the dissociation of no less than 0.5 ML 
of water molecules closes the band gap, resulting in a metallic interface. Figure |6tb) shows the 
details of the electronic states of water molecules and the LAO surface layer near the Fermi level. 
Since the VBM of the LAO surface layer and the water molecules exceeds the STO CBM, a part of 
charge transfers from the surface into the STO conduction band. The carrier density is equal to the 
density of the empty states in the valence band shown in Figure [6tb). As shown in Figure |6tc), the 
carrier density increases with increasing the coverage of dissociated water molecules. The carrier 
density is L3 x 10^^ cm"^ for c(l, 3/4) , and reaches to the highest value, L9 x 10^^ cm"^, for 
c(l,l). 
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FIG. 6. (Color online) (a) Layer-resolved projected densities of states (LPDOSs) for the configurations of 1 
ML coverage. Grey-filled lines represent the LPDOSs for fully molecular adsorption c(l, 0) , green, blue, 
and red solid lines represent the LPDOSs for c(l, 1/2), c(l, 3/4), and c(l, 1). (b) Details of LPDOSs of 
water molecules and the LAO surface layer near the Fermi level for c(l, 1/2), c(l, 3/4), and c(l, 1). Lines 
are colored as in (a). Vertical solid lines denote Fermi levels. The VBM of STO bottom layer is taken to 
eV. (c) Carrier densities for the configurations of 1 ML coverage with different coverages of dissociated 
water molecules (6). 



VI. DISCUSSIONS 



The above calculated results are schematically illustrated in Fig.|7l The electronic properties at 
the LAO/STO interface are nearly not affected by molecular adsorption of water molecules at the 
LAO surface. Dissociative adsorption shifts up the p states of the O atoms (like Os) neighboring 
the hydroxyl Ow-H on the Al atom, and induces a state originated from the hydroxyl Ow-H at 
the VBM of the LAO surface. Once the coverage of dissociated water molecules reaches to 1/2 
ML the VBM of the surface exceeds the STO CBM, resulting in a charge transfer from the sur- 
face to the conduction band of the STO. In ambient atmosphere at room temperature c(l, 3/4) is 
ground state and has a lower carrier density, while c(l, 1) is metastable and has a higher carrier 
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FIG. 7. (Color online) Schematic atomic structures and band diagrams of the LAO/STO interface with 
water adsorption on the surface, (a), (b) Molecular adsorption, and (c), (d) dissociative adsorption. Shade 
areas are the valence bands of the STO and LAO. Red and pink filled curves present O p states of the LAO 
surface layer, blue filled curves represent the states of water molecules. 

density. This implies that a structural transition between these two stable adsorption configura- 
tions can induce a conductance change at the LAO/STO interface. This picture provides us a 
way to interpret the mechanism of the IMT observed by Cen et al. . But, before elucidating it we 
should first clarify so-called insulator state in their experiment. In fact, before the AFM tip with 
positive voltages scanning the surface the interface already presents a conductance about 6.9 nS, 
which was referred to the insulating state by Cen et al, and after the AFM tip scanning the surface 
the conductance rises to about 7.6 nS. The ratio of conductance change is about ten percent. In 
Comparison with about five orders of magnitude resistance difference between the insulating and 
metallic interfaces reported by Thiel et al~, the interface with 6.9 nS conductance should be rather 
regarded as metallic. Considering this point and associating the conductance change observed in 
Cen et aUs experiment with the difference of carrier densities of two stable adsorption config- 
urations, c(l, 3/4) and c(l, 1), we propose a mechanism as the following. Before the AFM tip 
scanning the LAO surface the configuration of water molecules lies at the thermodynamic ground 
state c(l, 3/4), the interface shows a relatively low conductance. After the AFM tip with positive 
voltages scanning the LAO surface all water molecules in the scanned area dissociate, i.e., the 
configuration of water molecules in the scanned area becomes c(l, 1), thus the interface shows a 
higher conductance. As for how the AFM tip dissociates water molecules, a recent study on water 
molecule dissociation at a MgO surface by means of STM tip scanning has provided a possible 
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answer—, in which it was demonstrated that a tunneUing current passing through the adsorbed wa- 
ter molecule can induce a lateral vibration model of water molecule, consequently facilitating the 
dissociation. Therefore, it is highly possible that in Cen et aUs experiment water molecules in the 
area scanned by the AFM tip were fully dissociated due to a tunnelling current passing through 
them. 

Moreover, according to our calculations the carrier density change at the LAO/STO interface 
induced by different water adsorptions can only happen in the case of three LAO layers, being in 
good agreement with Cen et aUs experimental result. For the interface with two LAO layers the 
calculated gap between the VBM of the LAO surface layer and the STO CBM is 1 . 1 eV, which 
is too large to be closed by water molecule adsorption because the VBM up-shift due to water 
adsorption is only about 0.6 eV. In the case of four LAO layers, the interface is already metallic 
and has a carrier density of 2 x 10^^ cm"^ as reported by Thiel et al^, which may smooth over the 
carrier density change due to the different adsorption configurations. 

Noticeably, besides the aforementioned mechanism of conductance change at the LAO/STO 
interface, our results also imply another scheme to realize a sharper switch between insulating 
and metal states at the LAO/STO interface with three LAO layers by adsorbing and desorbing 
water molecules at the LAO surface. In this scheme the insulating state corresponds to the case 
of bare surface and the metallic state corresponds to an adsorption configuration of 1 ML of water 
molecules adsorbed at the surface. The insulator-metal transition can be realized by exposing the 
bare surface to a water vapor and the reverse transition can be realized by desorbing the water 
molecules from the surface. 

VII. SUMMARY 

In summary, we have studied electronic and structural characteristics of water molecule ad- 
sorption on the LAO surface of the n-type LAO/STO interface by using first-principles electronic 
structure calculations. We found that in an ambient atmosphere at room temperature the adsorption 
configuration of 1 ML of water molecules including 3/4 ML of dissociated molecules is most sta- 
ble and the configuration of 1 ML of dissociated water molecules is metastable. For the LAO/STO 
interface with three LAO layers, dissociation of no less than 0.5 ML of adsorbed water molecules 
results in a metallic interface, and the carrier density increases with increasing the coverage of dis- 
sociated water molecules. A structural transition between the ground-state and a metastable water 
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adsorption configurations can induce a carrier density change at the LAO/STO interface. Our find- 
ings provide new insights into creating and modulating electron carriers at the LAO/STO interface, 
which hopefully prompt the application of oxide interface electronics. In addition, our results also 
suggest that water molecule adsorption from the ambient atmosphere should be paid particular at- 
tention for experimentally studying electronic properties of metal oxides, especially for the metal 
oxides with a small band gap it may dramatically change some fundamental electronic properties 
of the materials, for instance, insulator-metal transition. 
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